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Abstract 

A search for events containing two or more high-transverse-momentum isolated 
leptons has been performed in ep colhsions with the ZEUS detector at HERA 
using the full collected data sample, corresponding to an integrated luminosity of 
480 pb~^. The number of observed events has been compared with the prediction 
from the Standard Model, searching for possible deviations, especially for multi- 
lepton events with invariant mass larger than 100 GeV. Good agreement with 
the Standard Model has been observed. Total and differential cross sections for 
di-lepton production have been measured in a restricted phase space dominated 
by photon-photon collisions. 
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1 Introduction 



The production of multi-lepton final states in electron-proton coUisionq^ is predicted 
within the framework of the Standard Model (SM). At HERA energies, the production 
cross sections are small for high transverse momenta, px, of the produced leptons and, 
along with the distributions of the kinematic quantities, can be calculated with high accu- 
racy in the SM. Therefore contributions from beyond the SM could either be observed as 
an increase of the visible cross sections or as a deviation from the predicted distributions. 

Multi-lepton final states were searched for by the HI Collaboration [1] using a luminosity 
of 463 pb~^. The observed overall numbers of di- and tri-lepton events were in good 
agreement with the SM predictions. However, some events with large transverse momenta 
were observed, exceeding SM predictions in this region. 

The analysis presented here is based on a luminosity of 480 pb~^ collected by the ZEUS 
experiment. Events with two or more high-p^ leptons (electrons or muons) were searched 
for and the total yields and distributions of kinematic variables were compared to SM 
predictions. In addition, the total visible and differential cross sections for di-lepton 
production were measured in the photoproduction regime, in which the incoming electron 
has small squared momentum transfer, < 1 GeV^. 



2 Experimental set-up 

The analysed data were collected between 1996 and 2007 at the electron-proton col- 
lider HERA using the ZEUS detector. During this period HERA operated with an elec- 
tron beam energy of 27.5 GeV and a proton beam energy of 820 GeV and, from 1998, of 
920 GeV, corresponding to centre-of-mass energies of 300 GeV and 318 GeV, respectively. 

A detailed description of the ZEUS detector can be found elsewhere [2]. A brief outline 
of the components that are most relevant for this analysis is given below. 

Charged particles were tracked in the central tracking detector (CTD) [3], which operated 
in a magnetic field of 1.43 T provided by a thin superconducting solenoid and covered the 
polar- angl^ region 15° < 6 < 164°. Before the 2003-2007 running period, the ZEUS 

^ Here and in the following, the term "electron" denotes gencrically both the electron (e~) and the 
positron (e"*"). 

^ The ZEUS coordinate system is a right-handed Cartesian system, with the Z axis pointing in the 
proton beam direction, the Y axis pointing up and the X axis pointing towards the centre of HERA. 
The polar angle, 0, is measured with respect to the proton beam direction. The coordinate origin is 
at the nominal interaction point. 
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tracking system was upgraded with a silicon microvertex detector (MVD) [4]. The high- 
resolution uranium-scintillator calorimeter (CAL) [5] consisted of three parts: the forward 
(FCAL), the barrel (BCAL) and the rear (RCAL) calorimeters. The smallest subdivision 
of the CAL was called a cell. The muon system consisted of rear, barrel (R/BMUON) [6] 
and forward (FMUON) [2] tracking detectors. The B/RMUON consisted of limited- 
streamer (LS) tube chambers placed behind the BCAL (RCAL), inside and outside a 
magnetised iron yoke surrounding the CAL. The barrel and rear muon chambers covered 
polar angles from 34° to 135° and from 135° to 171°, respectively. The FMUON consisted 
of six trigger planes of LS tubes and four planes of drift chambers covering the angular 
region from 5° to 32°. The muon system exploited the magnetic field of the iron yoke 
and, in the forward direction, of two iron toroids magnetised to ~ 1.6 T to provide a 
measurement of the muon momentum. 

The luminosity was measured using the Bethe-Heitler reaction ep e'jp by a luminos- 
ity detector which consisted of a lead-scintillator [7-9] calorimeter and, in the 2003-07 
running period, an independent magnetic spectrometer [10]. The fractional systematic 
uncertainty on the measured luminosity was 2.5%. 

The integrated luminosity of the samples corresponds to 480 pb"^ for events in which a 
search of electrons but no muons was carried out (electron channel) and to 444 pb~^ for 
events in which a search for either muons or electrons was carried out (muon channel). The 
slight difference in integrated luminosity is due to the requirement of a good performance 
of the detector components involved in the search. 

3 Standard Model processes and Monte Carlo simu- 
lation 

To evaluate the detector acceptance and to provide simulations of signal and background 
distributions, Monte Carlo (MC) samples of signal and background events were generated. 

The SM predicts that isolated multi-lepton final states are predominantly produced by 
two-photon interactions, 77 — > /"''/". The Grape MC event generator [11] was used to 
simulate these processes. It also includes contributions from 7Z and ZZ interactions, 
photon internal conversions and virtual and real Z production. It is based on the elec- 
troweak matrix elements at tree level. At the proton vertex, three contributions were 
considered: elastic, where the proton stays intact; quasi-elastic, where a resonant state 
is formed; and inelastic, where the proton interacts via its quark constituents. At the 
electron vertex, all values of were generated, from — GeV^ (photoproduction) 
to the deep inelastic scattering (DIS) regime. The uncertainty on the Grape predictions 
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was taken to be 3% [1]. 

The Drell-Yan process from resolved photon events, in which the photon fluctuates into 
a qq pair, and the lepton pair is produced from the interaction between a quark in the 
proton and one of the quarks from the photon, is not included. However this is expected 
to be negligible in the investigated kinematic regime [12]. 

The dominant SM background to topologies in which at least one electron is identified 
comes from neutral current (NC) DIS and QED Compton (QEDC) events. In NC (ep — 
eX) events, the scattered electron is identified as one of the electrons of the pair and 
hadrons or photons in the hadronic system X are misidentified as a further electron. In 
QEDC events (ep — > e7X), the final-state photon may convert into an e^e~ pair in the 
detector material in front of the CTD and typically one of these two electrons is identified 
as the second electron of the pair. 

The NC DIS and QEDC events were simulated with the Djangoh [13] and Grape- 
CoMPTON [11] MC programs, respectively. The absolute predictions of the Grape- 
COMPTON MC were scaled by a factor 1.13 in order to correct imperfections in the 
simulation of the dead material between the beampipe and the CTD. The uncertainty on 
this factor was taken as a source of systematic uncertainty. For the final-state topologies in 
which an electron and a muon were found, the background from SM di-tau pair production 
was estimated using the Grape MC program. 

Standard Model processes such as vector-meson (T, charmonium) and open heavy-flavour 
(charm and beauty) production were studied using the DiFFVM [14] and Pythia [15] MC 
programs and were found to be negligible. 

The generated events were passed through a full simulation of the ZEUS detector based 
on the GEANT [16] program versions 3.13 (1996-2000) and 3.21 (2003-07). They were 
then subjected to the same trigger requirements and processed by the same reconstruction 
program as the data. 

4 Event selection 
4.1 Online selection 

Events with two or more leptons in the final state were selected using the ZEUS three-level 
trigger system [2,17,18]. 

To select electrons, a significant energy deposit was required in the electromagnetic 
calorimeter and at least one good track in the central detectors had to be present. In 
addition, two other trigger chains were used: the first, dedicated to NC DIS selection. 
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requiring the detection of an electron with an energy E'^ > 4 GeV; the second, dedi- 
cated to the selection of events with high transverse energy deposited in the calorimeter 
{Et > 25 GeV). 

To select muons [19], a candidate was identified as a central track measured in the CTD 
matched to an energy deposit in the CAL and to a segment in the barrel or rear inner 
muon chambers. 

4.2 Electron identification 

The following criteria were imposed to select electrons in the offline analysis: 

• electron identification — an algorithm [20] which combined information from the en- 
ergy deposits in the calorimeter and, when available, tracks measured in the central 
tracking detectors was used to identify the electron candidates. Electron candidates in 
the central region (20° < 9'^ < 150°) were required to have energy greater than 10 GeV 
and a track matched with the energy deposit in the calorimeter. The matched track 
was required to be fitted to the primary vertex and to have a momentum of at least 
3 GeV and a distance of closest approach between the energy deposition and the track 
of less than 8 cm. Forward electrons (5° < 6*^ < 20°) were also required to have an 
energy greater than 10 GeV while, for electrons in the rear region (150° < 6*^ < 175°), 
the energy requirement was decreased to 5 GeV; 

• isolation — to ensure high purity, each electron candidate was required to be isolated 
such that the total energy not associated with the electron in an 77 — </> cone of radius 0.8 
centred on the electron was less than 0.3 GeV. This requirement was complemented, 
for electrons in the central region, by the request that no track with p-r > 1 GeV, other 
than the matching track, was contained in an 77 — cone of radius 0.4 centred on the 
electron; 

• QEDC background reduction — for the data collected in 2003-07, each track associ- 
ated with an electron candidate was required to have at least two hits in the MVD. 
This requirement removed photon conversions in the material between the MVD and 
the CTD. 

4.3 Muon identification 

The following criteria were imposed to select muons in the offline analysis: 

• muon identification — at least one muon candidate in the event was required to be 
reconstructed by the rear, barrel or forward muon chambers, matched to a track and to 
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an energy deposit in the calorimeter. In the case when only one muon in the event was 
reconstructed by the muon chambers, additional muons were also selected with looser 
criteria, by requiring a track pointing towards a calorimeter energy deposit compatible 
with that from a minimum ionising particle (mip). Each muon candidate was required 
to be associated with a track fitted to the primary vertex. The muon momentum was 
reconstructed using the central tracking devices, complemented with the information 
from the FMUON when available. The muon was required to have pi^ > 2 GeV, and 
to lie in the angular region 20° < < 26° (FMUON), 35° < < 160° (B/RMUON), 
20° < < 160° (mip); 

• isolation — to ensure high purity, each identified muon was required to be isolated 
such that only the matching track was contained in an 77 — cone of radius 1.0 centred 
on the muon. This cut, harder than in the electron selection, was used to reject 
background events in which a muon was found very close to a hadronic system, in 
particular in the e/x channel; 

• cosmic-muon background reduction — the reconstructed primary vertex had to be con- 
sistent with the HERA beam-spot position. If two muons were found, the acollinear- 
ity angle, Q, between the two muons had to satisfy cosQ > —0.995. For events with 
cosQ < —0.990, additional CAL timing cuts were applied. 

4.4 Event selection and classification 

The final event selection required the event vertex to be reconstructed with |Zvtx| < 
30 cm. At least two leptons, electrons or muons, had to be reconstructed in the central 
part of the detector (20° < 6*' < 150°). One of the leptons had to have > 10 GeV and 
the other > 5 GeV. Additional leptons identified as described in Sections 14.21 and 14.31 
could be present in the event. No explicit requirement on the charge of the leptons was 
imposed. According to the number and the flavour of the lepton candidates, the events 
were classified into mutually exclusive samples. 

For the measurement of the production cross section of e^e~ and fi~^fi~ pairs in the 
photoproduction regime, the cut {E — Pz) < 45 GeV was applied. This quantity was 
reconstructed in the electron case as 



where the sum runs over the corrected energies, E'i°^^ , of the CAL clusters and, in the 
muon case, as 



E-Pz = Y^ E,{1 - COs(^i)) - J2 ^mip(l - COs(^„iip)) + J2 ^muon(l " COs(^„)), (2) 




(1) 



mip 



muon 
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where Ei is the energy of the i^^ CAL cell and the [E—Pz) of the CAL mip was replaced by 
that of the muon track. This requirement selects events in which the scattered electron 
was lost in the beampipe and corresponds to a cut of < 1 GeV^ and on the event 
inelasticity, y = {E — Pz)/2Ee < 0.82, where E^ is the electron beam energy. The 
background from NC DIS and QEDC events is negligible in this sample, which will be 
referred to as the 77 sample in the following. 



5 Systematic uncertainties 



The following sources of systematic uncertainties were considered; the effect on the total 
visible cross section is given: 

• the muon acceptance, including the B/RMUON trigger, the reconstruction and the 
muon identification efficiencies, is known to about 7% from a study based on an 
independent elastic di-muon sample [21], resulting in an uncertainty of for 
muons; 

• the uncertainty on the efficiency of the CTD part of the trigger chain was esti- 
mated from a study based on an independent sample of low-multiplicity low-Q^ DIS 
events [22], resulting in an uncertainty of +5% for electrons and ±5% for muons; 

• the CAL energy scale was varied by its uncertainty of 3%, resulting in an uncertainty 
of (^3%) for electrons and negligible for muons; 

• the uncertainty on the efficiency of the CAL part of the muon trigger (±3%) and of 
the mip finder (±2%) resulted in an uncertainty of ±4% for muons; 

• the uncertainty on the measurement of the hadronic system was evaluated by using an 
alternative reconstruction oi E — Pz, resulting in an uncertainty of — L8% for electrons; 

• the scaling factor of the QEDC MC was varied between 0.95 and 1.31, as allowed by 
the comparison with a QEDC-enriched data sample, resulting in a negligible effect for 
both electrons and muons. 

The total systematic uncertainty was obtained by adding the individual contributions 
in quadrature. A 2.5% overall normalisation uncertainty associated with the luminosity 
measurement was included only in the systematic uncertainty of the total visible cross 
section. 
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6 Results 



The number of selected events in the data are compared to SM predictions in Table [H 
The following different di- and tri-lepton topologies are listed: ee, fifi, e/i, eee and e/i/i. 
The observed number of events is in good agreement with the predictions of the SM, 
according to which the NC DIS and QEDC processes give a sizeable contribution to the 
ee channel. Most of the events contributing to the e/i topology are predicted to come from 
di-muon production at high Q^, in which the beam electron is scattered at large angles 
and is therefore seen in the detector, while one of the muons is outside the acceptance 
region. A small contribution to this channel (~ 2 events) is predicted to come from di-r 
production, while the NC DIS background constitutes ~ 10% of the sample. 

Three four-lepton events, 2 in the ee/x/z and 1 in the eeee channel, were observed, to be 
compared to a SM expectation of ~ 1. The contributions from true four-lepton events 
are not included in the SM predictions and are expected to be small. Events with other 
multi-lepton topologies were searched for, but none was found. 

In Tables [T] and [2] the NC DIS and QEDC background contributions are given as limits 
at 95% confidence level (C.L.) when none or few events were selected from the back- 
ground MC samples. When this is done, the central value of the total SM prediction is 
determined as the most probable value (mode) of the convolution of the Gaussian signal 
distribution with the poissonian background distributions, and the uncertainty on the 
total SM prediction is determined by taking the 68% C.L. interval. 

Two events, one with three electrons in the final state and one with two muons and an 
electron, passing the three-lepton selection cuts, are shown in Fig. [H 

6.1 Kinematic distributions 

The distributions of the mass of the two highest-pT leptons in the event, M12, and of the 
scalar sum of the transverse momenta of all the identified leptons in the event, YIPt^ ^'^^ 
shown in Figs. [2] and [3] for all the observed di- and tri-lepton topologies, and are compared 
to SM predictions. The SM gives a good description of the data. In the mass region 
between 80 and 100 GeV, which is sensitive to production, 7 events were observed in 
the data, compatible with the predictions from the SM of ~ 9 events, including ~ 1 event 
from real production. 

The high-mass and high-^p^ regions are particularly sensitive to possible contributions 
from physics beyond the SM. The event yields for M12 > 100 GeV for all the observed 
di- and tri-lepton channels are summarised in Table O In the electron channels, 3 events 
at high masses are observed, to be compared with a SM prediction of 2.5. Two of these 
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events are observed in the eee topology, for which the SM expectation is 0.7. No event 
with Mi2 > 100 GeV is seen in the muon channels. The event yield for XIPt > GeV, 
combined for all the lepton topologies, is summarised in Tabled Two events at high-^p^^ 
are observed, to be compared with a SM prediction of ~ 1.6. 

The distributions of M12 and YIPt^ combined for all the di- and tri-lepton topologies, are 
shown in Fig. HI Also in this case, the data are well described by the SM predictions. 

6.2 Cross sections 

Total visible and differential cross sections for di-electron and di-muon production were 
determined in the kinematic region defined by: 

p!^ > 10 GeV, pf^ > 5 GeV, 20° < O^^'^ < 150°, < 1 GeV^ y < 0.82. 

The cross sections are given at ^/s = 318 GeV: the small (~ 5%) correction needed for 
the 1996-97 data sample was extracted from the MC. The effect of final-state radiation 
on the cross section was checked and found to be negligible. 

The total visible cross sections, corrected for acceptance, were measured to be 

a(77 ^ e+e") = 0.64 ± 0.051°;°^ pb (3) 

for the electron channel, and 

a(77^/iV) = 0.58 ± 0.07^^6 pb (4) 

for the muon channel. 

Since the muon and electron cross sections differ only marginally, they were combined 
in a single measurement, evaluated as the weighted mean of the two [23], assuming the 
systematic uncertainties to be uncorrelated. The systematic uncertainties of each mea- 
surement were symmetrised before the combination, by taking as systematic uncertainty 
the largest between the negative and the positive. The total visible cross sections are 
shown in Table HI compared with the SM predictions. 

Differential cross sections as a function of the invariant mass, M12, the transverse mo- 
mentum of the highest-pT lepton, p^, and the scalar sum of the transverse momentum 
of the two leptons, "^Px, are shown in Fig. [5], separately for electrons and muons. The 
di-electron, di-muon and combined cross sections are summarised in Table |5l The com- 
bination was done as described for the total visible cross section. Good agreement is 
observed between the data and the SM predictions. 
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7 Conclusions 



Events with two or more isolated leptons with high transverse momentum were observed 
using the full data sample taken with the ZEUS detector at HERA. The total number 
of multi-lepton events for different lepton configurations as well as their px and mass 
distributions were studied. No significant deviations from the predictions of the SM were 
observed. In addition, the total visible and differential cross sections for the e~^e~ and 
^~ signatures were measured in photoproduction and were observed to be in good 
agreement with the SM predictions. 
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ZEUS (£ = 480 pb"^) 



Topology 


Data 


Total SM 


Multi-lepton Production 


NC DIS 


Compton 


ee 


545 


563^29 


429^21 


74 ±5 


60± 10 




93 


106±12 


106±12 


< 0.5 




ejJL 


46 


42±4 


371^ 


4.5±1.2 




eee 


73 


751^ 


731^ 


< 1 


< 3 


e/i/i 


47 


48±5 


48±5 


< 0.5 




eeee 


1 




0.6±0.1 


< 0.4 


< 1 


ee/i/i 


2 


0.5^°:? 


0.4±0.1 


< 0.5 




All 4 leptons 


3 


1 4+0.7 


1.0 ±0.2 


< 1.4 


ee (77 sample) 


166 


185+1^ 




1.4±1.0 


1.4 ±0.6 


/i/i (77 sample) 


72 


°^-10 


9 

0^-10 


< 0.5 





Table 1: The observed and predicted multi-lepton event yields for the ee, nn, 
e/i, eee, e/x/i, eeee and eefifj, event topologies; the event yields for the ee and fifj, 
topologies in the 77 samples. The quoted uncertainties consist of model uncertain- 
ties, MC statistical uncertainties and systematic experimental uncertainties added 
in quadrature. Limits at 95% C.L. are given when none or few events were selected 
from the background MC samples. The central value and the uncertainty on the 
total SM predictions are in these cases determined as explained in the text. 



ZEUS {C = 480 pb~^) 



Topology, 
M12 > 100 GeV 


Data 


Total SM 


Multi-lepton Production 


NC DIS 


Compton 


ee 


1 


1.7 ±0.2 


0.9 ±0.1 


0.2 ±0.1 


0.6 ±0.1 


/i/i 





0.4 ±0.1 


0.4 ±0.1 


< 0.01 




efj, 





'-'•"-'"-0.01 


0.05 ±0.02 


< 0.02 




eee 


2 


0.7 ±0.1 


0.7 ±0.1 


< 0.01 


< 0.02 


ejjLjjL 





0.18 ±0.05 


0.18 ±0.05 


< 0.01 





Table 2: The observed and predicted high-mass, M12 > 100 GeV , multi-lepton 
event yields. The invariant mass was calculated using the two highest-px leptons. 
The quoted uncertainties consist of model uncertainties, MC statistical uncertain- 
ties and systematic experimental uncertainties added in quadrature. Limits at 95% 
C.L. are given when none or few events were selected from the background MC 
samples. The central value and the uncertainty on the total SM predictions are in 
these cases determined as explained in the text. 
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ZEUS (£ = 480 pb"^) 



J2Pt> 100 GeV 


Data 


Total SM 


Multi-lepton Production 


NC DIS 


Compton 


2 


1.56 ±0.15 


1.16 ±0.13 


0.05 ±0.02 


0.35 ±0.06 



Table 3: The observed and predicted high-Y^p^j^ multi-lepton event yields for all 
topologies combined, where X^Py was calculated using all the leptons in the event. 
The quoted uncertainties consist of model uncertainties, MC statistical uncertain- 
ties and systematic experimental uncertainties added in quadrature. 



ZEUS (£ = 480 pb~^) 



77 sample 


96-07 (^■\^\ 
^DATA (pb) 


c^SM (pb) 


ee 

Combined 


0.64 ± Q.^htm 
0.59 ± 
0.63 ±0.04 ±0.03 


0.71 ±0.02 
0.69 ±0.02 
0.70 ±0.02 



Table 4: Total cross section for 77 ee and 77 — /x/i samples, combined as 
explained in the text, compared with the predictions from the Grape MC. 
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ZEUS (£ = 480 pb"^) 



Bin ( GeV) 


<TA(fb/GeV) 


0"SM 

(fb/GeV) 


e+e 




Combined 


15 < Mi2 < 25 
25 < Mi2 < 40 
40 < Mi2 < 60 
60 < Mi2 <100 


23.0 ±3.2 t^:^ 

19.1 ±2.0 1^:^ 
3.8 ±0.8 toi 

0.15 ± 0.11 


32.7 ±5.0 til 
12.6 ±2.6 
2.5 ±1.0 t'ol 
0.21 ±0.21 


25.4 ± 2.7 ± 1.5 
16.3 ± 1.6 ± 1.1 
3.3 ±0.6 ±0.3 
0.17 ±0.10 ±0.03 


30.4 ± 1.0 
19.8 ±0.7 
3.0 ±0.1 
0.26 ±0.02 


10 <p^ < 15 
15 <p^ < 20 
20 < < 25 
25 <p^^ < 50 


90.7 ±8.4+^;^ 
26.4 ±4.3+?:^ 

4.2±1.7toi 
0.90 ± 0.37+[J:?^ 


94 ± 12+^^ 
10.7 ±4.4+};? 
8.9±4.0^J:^ 
0.70 ± 0.50t[J:J^^ 


91.6 ±6.9 ±5.3 
18.1 ±3.1 ± 1.3 
5.0 ± 1.6 ±0.7 
0.82 ±0.29 ±0.08 


103.2 ±3.3 
23.7 ±0.9 
7.3 ±0.4 

0.88 ±0.06 


15 < J2Pt < 25 
25 < J2Pt < 40 
40 < J2PT < 60 
60 < J^Pt < 100 


36.7 ± 4.01?:^ 
15.8± 1.9+ii 
1.24±0.44+°:2i 
0.092 ± 0.092+°;°^! 


38.4±5.5l^:^ 
11.2 ±2.61};? 
1.32 ± 0.76l|]:}| 
0.18±0.18t[!:J^^ 


37.2 ±3.2 ±2.2 
14.0 ± 1.5 ±0.9 
1.26 ±0.38 ±0.16 
0.11 ±0.08 ±0.02 


43.9 ± 1.4 
14.6 ±0.5 
1.69 ±0.11 
0.16 ±0.02 



Table 5: Differential cross section as a function of the invariant mass, M^, the 
transverse momentum of the highest-pT lepton, p^^, and the scalar sum of the trans- 
verse momentum of the two leptons, XIPt' /^'^ di-lepton events in the kinematic 
region defined in the text, compared with the predictions from the Grape Monte 
Carlo. The results are shown separately for the ee and fifi samples, as well as for 
the combined sample. 
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Figure 1: a) An event with three electron candidates in the ZEUS detector. The 
invariant mass of the two highest-px electrons is M12 = 113 GeV ; the correspond- 
ing transverse momenta are given above, b) An event with two muons and an 
electron candidate in the ZEUS detector. The invariant mass of the di-muon pair 
is 77.5 GeV ; the corresponding transverse momenta are given above. 
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Figure 2: Distributions of the invariant mass of the two highest-px leptons for 
the different multi-lepton topologies: ee, eee, /x/i, e/i/i, e/x. The ZEUS data are 
displayed as the full dots. The errors on the data are given by the square root of the 
number of events in each bin. The SM predictions are represented as the solid line 
and are obtained by summing the contributions of di-lepton production, NC DIS, 
QED Compton events, and, for the en channel, di-tau production. The error band 
represents the systematic uncertainty on the SM predictions. 
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Figure 3: Distributions of the sum of the transverse momenta of the leptons for 
the different multi-lepton topologies: ee, eee, fifi, cfifi, cfi. The ZEUS data are 
displayed as the full dots. The errors on the data are given by the square root of the 
number of events in each bin. The SM predictions are represented as the solid line 
and are obtained by summing the contributions of di-lepton production, NC DIS, 
QED Compton events, and, for the e/i channel, di-tau production. The error band 
represents the systematic uncertainty on the SM predictions. 
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ZEUS 




Figure 4: Distributions of (a) the invariant mass of the two highest-px leptons and 
(h) the sum of the transverse momenta of the leptons for all the individual lepton 
topologies combined: ee, eee, /i/x, e/x/i, e/i. The ZEUS data are displayed as the full 
dots. The errors on the data are given by the square root of the number of events 
in each bin. The SM predictions are represented as the solid line and are obtained 
by summing the contributions of di-lepton production, NC DIS, QED Compton 
events, and, for the efi channel, di-tau production. The error band represents the 
systematic uncertainty on the SM predictions. 
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Figure 5: Differential cross sections as a function of (a) the invariant mass of 
the lepton pair, M12, (h) the transverse momentum pj. of the highest-px lepton, 
and (c) the scalar sum of the transverse momenta of the two highest-px leptons, 
'Y^Px- The di-muon cross sections are shown as the open dots, while the full dots 
are the di- electron measurements, which have been displaced for clarity. The data 
are compared with the predictions of the Grape Monte Carlo. The full error bars 
are the quadratic sum of the statistical (inner part) and systematic uncertainties. 
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